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Abstract 

We predict a condensation phenomenon in an overall neutral system, consist- 
ing of a single charged plate and its oppositely charged counterions. Based 
on the "two-fluid" model, in which the counterions are divided into a "free" 
and a "condensed" fraction, we argue that for high surface charge, fluctua- 
tions can lead to a phase transition in which a large fraction of counterions 
is condensed. Furthermore, we show that depending on the valence, the con- 
densation is either a first-order or a smooth transition. 
61.20.Qg, 61.25.Hq, 87.15.Da 
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I. INTRODUCTION 



Electrostatic interactions control the structure, phase behavior, and function of 
macroions in aqueous solutions |IJ. The macroions may be charged membranes, stiff poly- 
electrolytes such as DNA, or charged colloidal particles. The fundamental description of 
these charged systems has been the Poisson-Boltzmann (PB) theory. However, it ignores 
fluctuations and correlations, which are important for the cases of low temperatures, highly 
charged surfaces, or multivalent counterions. These fluctuation and correlation effects, which 
have been the focus of recent theoretical efforts, may drastically alter the mean-field picture 
of PB theory [@-|5|]. For example, one surprising effect [f| is the attraction between two 
highly charged macroions, as observed in experiments || and in simulations @]. In this 
paper, we argue that correlation effects may lead to condensation of counterions onto an 
oppositely charged plate, whose surface charge becomes effectively renormalized. In par- 
ticular, the counterion valence plays an interesting role: for Z > Z c ~ 1.62 for typical 
system parameters (see below), we find a first-order phase transition in which a large frac- 
tion of the counterions is condensed, while for Z < Z c the condensation proceeds smoothly, 
implying that monovalent and divalent counterions exhibit qualitatively distinct behavior. 
This is in contrast with more familiar theories of counterion condensation ||, e.g. Manning 
condensation for charged rods, where the effective charge is continuously modified by the 
valence. 

Recall that for a single plate of charge density <r(x) = o~q5(z) immersed in an aqueous 
solution of dielectric constant e, containing point-like counterions of charge — Ze on both 
sides of the plate, PB theory predicts that the counterion density |] 

c(z) = (1) 

V ; 2tiZ 2 I b {\z\ + A) 2 V ; 

decays to zero algebraically with a characteristic length A = e/ (kIbZoq), where Ib = ~ 
7 A is the Bjerrum length in water at room temperature, ks is the Boltzmann constant, and 
T is the temperature. This Gouy-Chapman (GC) length A defines a sheath near the charged 
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surface within which most of the counterions are confined. Typically, it is on the order of few 
Angstroms for o~ ~ e/100 A" 2 . Note that Eq. ([!]) implies that at zero temperature all of the 
counterions would collapse onto the charged plane. However, for high surface charge (or low 
temperature) Z 2 Ib ^> A, fluctuation and correlation corrections can become so large that 
the solution Eq. ([I]) to the PB equation is no longer valid Therefore, we might expect a 
quantitative deviation from the conclusion above. Indeed, as pointed out by Netz et al. |§, 
a perturbative expansion about the PB solution breaks down in this regime, as indicated by 
an unphysical (negative) counterion density in the one-loop approximation. Motivated by 
these observations, we propose a two-fluid model in which the counterions are divided into a 
free and a condensate fraction. The free counterions have the usual 3D spatial distribution, 
while the condensed counterions are confined to the two-dimensional charged plane, with a 
mean (2D) density n c . We treat the fraction of 2D condensed counterions r = Zen c /cr as 
a variational parameter, which is determined self-consistently by minimizing the total free 
energy of the system. 

It may be useful to illustrate the essential physics first by a simple picture. In the spirit 
of the two-fluid model, the 2D condensed counterions partially neutralize the charged plate, 
effectively reducing the surface charge density from o"o to en# = o"o — Zen c , where n c is their 
surface (2D) density. The free counterions can be modeled as a 3D ideal gas confined to 
a slab of thickness \r = IKitIbZur). At the Debye-Hiickel level, the free energy per unit 
area for the condensed counterions f2d( n c) can be written as || 



Phd{n c ) = n c {ln[n c a 2 } - l} + i J jln 
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where (3~ l = ksT and A^ = l/(27r/BZ 2 n c ) is the 2D screening length. The first term in 
Eq. (0) is the entropy and the second term arises from the 2D fluctuations. Note that the 
latter term is logarithmically divergent, which may be regularized by a microscopic cut-off 
q c ~ 2n/a, yielding /3A/ 2 d(^ c ) — — ln(27rAi)/a). The free energy of the free counterions 
f3d(n c ) consists of the entropy of a confined 3D ideal gas and the fluctuation free energy. 
The latter term may be estimated by using the fluctuation contribution to the free energy 



density from the 3D Debye-Hiickel theory |T(J and multiplying it by the thickness of the 
slab X R : 



Pfad(n e ) « c\ R {ln[ca 3 ] - l} - X R , (3) 
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where c = n R /(ZXfi) average (3D) concentration of the free counterions and k 2 = AttZ 2 Ibc 
is the 3-D screening length. Note that the second term scales as ~ — A^ 2 . This simple 
picture to estimate fzd{n c ) contains all the qualitative physics fill , which follow from the 



more precise analysis presented below. The total free energy in the two-fluid model is 
/(t) = f2d{r) + / 3 d(r). Minimizing f(n c ) to find n c , we obtain 



1 + r " ln (4)" ln 



-iy(l-r) = 0, (4) 



1 - r)% 

where the three dimensionless parameters: the order parameter r = Zen c /<jQ, the coupling 
constant g = Z 2 Ib/X, (where A is the bare GC length), and the reduced temperature 9 = 
a/{Z 2 lB)-, completely determine the equilibrium state of the system. It is straightforward 
to obtain the asymptotic solutions of Eq. (f|) corresponding to the uncondensed, T\ <C 1, 
and condensed, r 2 ~ 1, state of the counterions: T\ ~ g9exp 1 — and r 2 ~ 1 — 
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[7rexp(l)p f 2 -) 2 . For weak couplings, g <C 1, T\ is the only consistent solution. On the 
other hand, for large coupling g ^> 1, where fluctuation free energies dominate the system, 
n and r 2 are both consistent solutions for small 9, and a first-order transition takes place 
when fiji) = /(t 2 ). Thus, a large fraction of counterions is condensed if g exceeds some 
threshold value g > g Q . For an estimate, taking 9 = 0.02 (divalent counterions at room 
temperature) we find go ~ 1.757, corresponding to a surface charge of <7o ~ e/10 nm -2 . 

We emphasize that although there is a close analogy between our approach and the more 
familiar theory of counterion condensation, e.g. Manning condensation ||, the counterion 
condensation in our model has a different physical origin arising from charge fluctuations. 
In Manning condensation, the competition between entropy and electrostatics leads to an 
electrostatic potential at large distances that is independent of the charge density of the 
rod above the Manning threshold ||. In this sense, for the geometry of a charge plate, 



counterions are always "Manning condensed" at the PB level ||12|| . On the other hand, in 
our model, we take one step further by showing that when correlation effects are taken into 
account, a finite fraction of the counterions is condensed to form a 2D Coulomb gas onto 
the charged plate. This paper is organized as follows: In Sec. [TI]> we present in detail the 
two-fluid model and construct the total free energy of the system. In Sec. fTJ we present 
the central results of this paper, followed by an extensive discussion. 

II. COUNTERION FREE ENERGY IN THE "TWO-FLUID" MODEL 

To study the condensation more rigorously, we compute total free energy by mapping the 
problem into a field theory. Consider an overall neutral system consisting of counterions and 
an oppositely charged surface immersed in an aqueous solution. The surface charge density 
on the plate is ctq = en . We model the aqueous solution with a uniform dielectric constant e. 
This simplification allows us to study fluctuation and correlation effects analytically. In the 
spirit of the "two-fluid" model, we divide the counterions into a "condensed" and a "free" 
fraction. The condensed counterions are allowed to move only on the charged surface, while 
the free counterions distribute in the space on both sides of the plate. The electrostatic free 
energy for the whole system may be written as 

(3F e i = J d 2 rn c {r) {in [n c (r)a 2 ] - l} + J d 3 xp(x) {in [p(x)a 3 ] - l} 

+ ^ r d3x r dSx ,n c (r)5(z)n c (r>)5(z>) ^ r t ■ . p(x) P (xQ 
2 J J |x-x'| 2 J J |x-x'| 
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where a is the molecular size of the counterions, Ib = e 2 / (eksT) is the Bjerrum length, Z is 
the valence of the counterions, r is the in-plane position vector, and x = (r, z) . The first two 
terms in Eq. (|5]) are the two-dimensional entropy for the condensate and three-dimensional 
entropy for the "free" counterions, respectively, and the other terms represent the electro- 
static interactions of counterions in the system. In Eq. (||), the two-dimensional density of 



the condensed counterions is denoted by n c (r), the "free" counterions with 3D density by 
p(x), and the external fixed charges arising from the surface by n./(x) = tiq5{z). Within 
the Gaussian fluctuation approximation, we consider the spatial dependent fluctuations of 
the 2D density of condensed counterions about a uniform mean: n c (r) = n c + 5n c (r), and 
expand Eq. 



to second order in 5n c (r): 



f3F ei = n c {in 
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where ^4 is the area of the plane, J(x) = -Zp(x') — njj(x), and n^(x) = n/(x) — Zn c 5(z). 
Note that J(x) is linearly coupled to 5n c (r) in the above equation. Summing over all the 
2D fluctuations of the condensed counterions, i.e. 



-me 



J V8n c (r)e- pF « 



we obtain two terms in the effective free energy: TC e = F 2d + 'H^d- The first term F 2d is the 
free energy associated with the condensed counterions which can be written as 



pF 2d = n c {ln[n c a 2 ] - l} A + - lndet K 2d - - lndet[-V 
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where K 2rf (x,y) 
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•Vx + -^S(z) 8(x — y) is the 2D Debye-Hiickel operator and \d = 
l/(2TrZ 2 l B n c ) is the Debye screening length in 2-D. The first term in Eq. (|j) is the entropy 
and the second term arises from the 2D charge-fluctuations. Note that this fluctuation term 
can be evaluated analytically 0, with the result quoted in Eq. (0): 



Pf2d(n, 



The second term Ti^d is the electrostatic free energy for the "free" counterions, taking 
into account of the presence of the fluctuating condensate; to within an additive constant, 
it may be written as 



PH 3d = J rf 3 xp(x) {in [p(x)a 3 ] - l} + i Jd 3 xJ d 3 x'p(x)G 2d (x,x>(x') 

d 3 x0(x)p(x), (7) 



where 0(x) = / d 3 x' G 2 d(x, x') n^(x') is the "renormalized" external field arising from 
the charged plate. From Eq. (0), we can see that the presence of the condensate modifies the 
electrostatics of the free counterions in two ways. First, the condensate partially neutralizes 
the charged surface, effectively reducing the surface charge density from eriQ to en^ = e(no — 
Zn c ). Second, their fluctuations renormalize the electrostatic interaction of the system; 
thus, instead of the usual Coulomb potential, the free counterions and the charged surfaces 
interact via the interaction G 2 d(x, x'), which is the inverse (the Green's function) of the 2D 
Debye-Hiickel operator K 2d [fH|: 



Ad 



G 2d (x,x') = 4tt/ b Z 2 5(x - x') 



(8) 



where the second term in the bracket takes the fluctuating 2D "condensate" into account. 
Hence, in the limit n c — > or Ad — > 00, G 2 d(x, x') reduces to the usual Coulomb interaction 
G (x,x') = 4vr/ B Z 2 /|x-x / |. 

After a Hubbard-Stratonovich transformation fl4| , the grand canonical partition function 
for the free counterions can be mapped onto a functional integral representation: 2^[0] = 
7V j Vijj e^'ft with the effective Hamiltonian p] 



J rf 3 x { 1 V(x)[-V 2 ]V(x) + [V(x)] 2 - «a e *M+*« 



(9) 



where V( x ) is the fluctuating field, k 2 = AtiIbZ 2 e^/a 3 , fi is the chemical potential, and 
Mq 2 = det K 2 rf is the normalization factor. The minimum of the effective Hamiltonian, 



given by 
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8tp(x) 



1p = lpQ 



0, defines the saddle-point equation for ^(x), which reads 



VV(x) + kV v(x) = Anl B Zn R 5(z) + — 5(z) <p(x) 

Ad 



(10) 



in terms of the mean-field potential <p(x) = — i"0o( x ) ~~ <K X )- The solution to Eq. ( [TO] ) 
is <p(x) = 2 In f 1 + , which satisfies the boundary conditions: i) <p(0) = and ii) 



dip 
dz 



= IixIrZur, with k = IhIbZurI \pl. Thus, at the mean-field level, the distribution 

2=0 



of the free counterions 

p (x) = n 2 e-^ x) /4nl B Z 2 



2txI b Z\\z\ + \ R )i 

has exactly the same form as the PB distribution Eq. (p]), but with a renormalized GC length 
\ R ; = \^2/k = 1/ '(tvIbZiir). To obtain the mean- field free energy of the free counterions 
Fo(nji), we note that it is related to the Gibbs potential r o [0] = S[ipo,<p] by a Legendre 
transformation: F (nn) = r o [0] + fi f d 3 x/) (x). Solving for the chemical potential fi from 
its definition: /i = In ( ^Ar ) an< ^ us ^ n S ^he m ean-field solution y?(x), we find 

^oMM=-ln^— j- T , (11) 

where ^4 is the area of the charged plane. Note that Fq{jir) has the form of an ideal 
gas entropy of a gas with concentration tir/ (ZXr) confined to a slab of thickness Xr, the 
renormalized GC length. 

Next, to capture correlation effects, we must also include the fluctuations of the free 
counterions, thereby treating the "free" and "condensed" counterions on the same level. 
To this end, we expand the action S[tp,(j)] about the saddle-point V-'o(x) to second order in 
A^(x) = V(x) - Vb(x) 

S[<f>, ij] = S[<f>, M + l I d 3 * J d 3 y A^(x) K 3rf (x, y) A^(y) + • • • , (12) 

where the differential operator 
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K 3d (x,y) 



5(x-y), (13) 



Ad v ' (M + A R ) 2 

is the second variation of the action S[i[),<f)]. Note that the linear term in A^(x) does 
not contribute to the expansion since Vo( x ) satisfies the saddle-point equation Eq. (|HID. 
Performing the Gaussian integrals in the functional integral, we obtain an expression for the 
change in the free energy due to fluctuations of the free counterions: 

(3AF 3d = hndetK 3d -hndetK 2d , (14) 



where the second term comes from the normalization factor j\fo- To evaluate (3AF^ explic- 
itly, we first differentiate it with respect to l B by making use of the identity <51ndetX = 
Tr X -1 5 X to obtain 

d\ R f ^ 2G M (x,x) dX D fja^fr, Ur „ N ^ Ur „ N1 S(z) 
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|rf 3 x [G M (x,x)-G M (x,x)] 
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(15) 
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where G 2d (x, x) = / ^ Ma* 
3D free counterions. It satisfies 

which can be solved to yield: 



and G3d(x, x') is the Green's function for the 
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where 7 = Xr/Xd = 2 t/ (1— r). Note that the first term in G 3( i(x, x) is just the Coulomb self- 



energy G (0) = / ^ 



<2 2 q 2ttI r Z 2 



which must be subtracted. Inserting G 2 d(x, x) and G 3 ^(x, x) 



into Eq. fllED , we obtain 

1 d(3AF 3d _ Jx(7) dX R J 2 ( 7 )/7 <9A d 

A dl B AttX 3 r dl B AttX% dl B 

where the functions Ii ,2(7) are given by 



+ (self-energy), 



(18) 



X 1 (7) = -ln(l + 7 ) + 3 
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Because X 12 (7) are independent of l B , we can integrate Eq. ( [15] ) back to obtain f3AF 3d ; 
thus, the total free energy per unit area for the free counterions is 



Pf3d[T) ~ Z \2ZXr) Z 8nX R 8ttX d Xr 



(19) 



Incidentally, in the limit of vanishing density of the condensed counterions, n c — > (or 
Xr — > A), Xl(0) = and we obtain the fluctuation correction to the mean-field PB free 
energy: Af P b = -k B T/(8y/3X 2 ). This result may be understood physically as follows. 
According to PB theory, the counterions are confined to a slab of thickness A, and thus may 
be considered as an ideal gas with a 3D concentration of c ~ no/A. This implies that the 
inverse of the 3D "screening" length is k s ~ \[cTb ~ 1/A. Using the 3D Debye-Hiickel free 
energy (per unit volume) /3Af ~ — k^, the correction to the mean-field PB free energy (per 
unit area) scales like (3Af p b ~ — A • A~ 3 ~ — A -2 , in agreement with Eq. (Q). Therefore, the 
precise calculation leading to Eq. (plj) justifies the use of the simple picture to illustrate the 
physics behind the counterion condensation presented in the Introduction. We note finally 
that Eq. (|T9|) also contains additional couplings among the fluctuations of the "condensed" 
and "free" counterions. 

III. RESULTS AND DISCUSSION 

The central results of this paper follow from the minimization of the total free energy 
f( T ) = hd{ T ) + fzd (t), obtained respectively in Eqs. (§) and fljjD , with respect to the order 
parameter r. Fig. [I] summarizes the behavior of r as a function of the coupling constant 
related to the surface charge density g = Z 2 Ib/X and the reduced temperature 9 = a/ [ZHb)- 
For weak coupling g 1, where fluctuation corrections are negligibly small, the counterions 
prefer to be free to gain entropy; there are almost no condensed counterions so that r w 0. 
This is not surprising since PB theory is a weak-coupling theory which becomes exact as 
g — > 0. However, for higher surface charge density, where correlation effects become more 
important, the behavior of r depends crucially on 6. In particular, for 9 < 9 C w 0.0378, r 
displays a finite jump at go(9), e.g. g = 1.695 at 9 = 0.02. [This corresponds to divalent 
counterions at room temperature with <jq ~ O.le nm~ 2 .] Thus, the system exhibits a first- 
order phase transition, in which a large fraction of counterions is condensed (about 80%). 
The physical mechanism leading to this counterion condensation is the additional binding 
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arising from 2D charge-fluctuations, which dominate the system at lower temperatures. 
However, for 9 > 9 C the behavior of r is completely different; in this regime, there is no 
phase transition and the condensation occurs smoothly. Thus, the condensation transition 
is similar to the liquid-gas transition, which has a line of first-order transitions terminating 
at the critical point where a second-order transition occurs. In our case, the critical point is 
found to be r c ~ 0.4, g c ~ 1.605, and 9 C ~ 0.0378. Furthermore, if one takes Ib ~ 10 A, i.e. 
room temperature, and a ~ 1 A, it follows from the definition of 9 that there is a critical 



value of counterion valence Z c = ya/(/_e^ c ) — 1.62, below which no first-order condensation 
transition is possible. Therefore, divalent counterions behave qualitatively differently from 
monovalent counterions. In fact, significant differences between mono- and divalent ions are 
observed in various biophysical processes. 

We stress that fluctuation effects are crucial for this counterion condensation transition 
to occur. In fact, it may be viewed as a surface analog of the bulk transition discussed 
by Fisher and Levin ||16|| . These authors predicted a phase separation, where a strongly 
correlated, dense phase coexists with a weakly correlated dilute phase in an ionic system 
dominated by Coulomb interactions and charge fluctuations. In our case, the surface breaks 
the translational symmetry and similar phase separation occurs in its vicinity. Indeed, using 
Eq. (|i~9|) in the limit of r — > 0, i.e. without assuming the existence of 2D condensate, 
the system shows a thermodynamic instability at g ~ 4.4. The inclusion of an additional 
degree of freedom, i.e., allowing the counterions to condense, can only lower the total free 
energy, suggesting a phase transition in which the condensate ("liquid") near the surface 
coexists with the more dilute, delocalized counterion ( "gas" ) distribution. Indeed, a recent 
simulation [17] clearly shows that at low temperature, most of the counterions reside on 



the surface, consistent with our two-fluid picture. However, our calculation based on the 
Gaussian fluctuation theory may break down for very large g > 10. In this regime, a 
complementary treatment is considered by Shklovskii in Ref. 0], in which the condensed 
counterions are assumed to form a 2D Wigner crystal. That theory also predicts a strongly 
reduced surface charge and an exponentially large renormalized Gouy-Chapmann length, 
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qualitatively similar to our results. In contrast, by treating the fluctuations of the condensed 
and free counterions on an equal footing, we are able to capture the onset of the condensation 
(at g ~ 2), which bridges between the regime where PB theory is appropriate, g — > 0, and 
the very strong coupling regime, g — > oo PJTTfl. 



In summary, we have presented a new mechanism by which the counterions become 
condensed so as to neutralize the surface charge of a macroion. It has been known experi- 
mentally that an effective surface charge, which is always lower than the actual charge, must 
be introduced in order to fit experimental data to the PB theory fl8| . Thus, our theory offers 



a possible scenario to account for this experimental fact. In addition, for the case of two 
highly-charged surfaces, the PB repulsions between them are greatly reduced due to strong 
condensation, and the dominant interaction will be the charge fluctuation attractions. Thus, 
this condensation picture may also be crucial to understanding the like-charged attraction 
fT9|| . Furthermore, there are some recent experimental |2(| and simulation indications 



that are consistent with the predicted condensation effect. The experiments |2t| were per- 
formed with a monolayer of cationic surfactant where surface density of the surfactant and 
the counterion/salt density are controlled with high accuracy. The experiments measured 
a rapid neutralization (about 90%) of the charged surfactant monolayer by increasing its 
surface density (by about 10%). In some cases, a discontinuous neutralization process is 
observed ||20|| . Also, recent extensive simulation studies of uniformly charged surfaces per- 
formed in Ref. |21] reports two interesting observations: First, when g 3> 1 there appears a 



coexistence between two distinct counterion density distributions: an exponentially decaying 
distribution near the immediate vicinity of the charged surface and an algebraic decaying 
distribution far away from the the surface. Note that the exponentially decaying distribu- 
tion might be associated with our condensed counterions, which we have assumed to be 
a delta-function distribution. Secondly, the specific heat of the simulated system shows a 
pronounced hump in the region 10 < g < 100, though no rigorous proof of the condensa- 
tion transition from simulations (and experiments) has been obtained so far. Indeed, there 
remains some fundamental issues to be addressed in the future, for example, the role of 
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excluded volumes, the discreteness of the surface charge and its mobility, and higher order 
(beyond Gaussian) corrections. A recent calculation and simulation shows that charge dis- 



creteness also induces charge localization [p^j . Therefore, it is possible that these neglected 
effects may smooth out the first-order transition. However, we believe that a rapid variation 
of the condensation with the surface charge, reflecting the predicted effect, should remain. 
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FIG. 1. The fraction of condensed counterions r = Zen c /ao as a function of g = Z 2 Ib/X for 
different values of 9 = cl/{Z 2 Ib)- At low surface charge g <C 1, the counterion distribution is well 
described by PB theory since However, at high surface charge, correlation effects leads to 

a large fraction of counterion condensed. The condensation is first-order for 9 < 6 C and smooth for 
9 > 9 C , where 9 C pa 0.0378. The solid line 9 = 0.02 corresponds to divalent counterions, where a 
finite jump occurs at go ~ 1.7 or ctq ~ 0.1 e nm~ 2 . 
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